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Abstract 

In this paper we study the phase diagram of the disordered Ising ferromagnet. Within the 
framework of the Gaussian variational approximation it is shown that in systems with a finite 
value of the disorder in dimensions D = 4 and L) < 4 the paramagnetic and ferromagnetic 
phases are separated by a spin-glass phase. The transition from paramagnetic to spin-glass 
state is continuous (second-order), while the transition between spin-glass and ferromagnetic 
states is discontinuous (first-order). It is also shown that within the considered approximation 
there is no replica symmetry breaking in the spin-glass phase. The validity of the Gaussian 
variational approximation for the present problem is discussed, and we provide a tentative 
physical interpretation of the results. 
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1 Introduction 



Phase transitions in the random temperature Ising ferromagnets have been intensively studied 
theoretically, numerically and experimentally during the last decades. The theoretical interest has 
mainly been focused on the critical behavior in the vicinity the paramagnetic-ferromagnetic phase 
transition point in weakly disordered systems Renormalization group considerations show 
that if the temperature is not too close to T^, the critical behavior is essentially controlled by the 
fixed point of the pure system (so that disorder produces only irrelevant corrections), while in the 
close vicinity of Tc the critical behavior turns out to be different from that of the pure system and 
is characterized by a new universal (independent of the disorder strength) fixed point. 

On these grounds it is widely believed that the critical behavior of the disordered system is 
universal, and the strength of the disorder only affects the size of the critical region near Tc (but 
not the critical behavior itself). In other words, the critical behavior of systems with a finite value 
of the disorder must be the same as that in the weakly disordered ones. Most of the numerical 
simulations (in particular for the two-dimensional systems) support this idea, see e.g. 0, although 
some of the numerical results seem to indicate that the critical behavior can be non-universal and 
characterized by critical exponents depending on the disorder strength 0. 

In this paper we consider the problem of the phase transitions in the disordered Ising ferromagnet 
from a somewhat different point of view. Instead of studying of the critical behavior, we propose 
first to address a simpler point concerning the nature of the phases in such systems. We stress 
that the usual assumption that the random temperature Ising ferromagnet can only be in the 
paramagnetic or in the ferromagnetic state is, to the least, questionable, a point first raised by Ma 
and Rudnick 0]. Before studying such details as the critical properties one should first clarify what 
kind of phases and what kind of transitions can exist in such a system. 

A general reason for asking such a question comes from the fact that the saddle-point equations 
which describe the local minima of the disordered Hamiltonian in the (supposed) paramagnetic 
region have an exponentially large number of solutions. Physically this situation is quite clear: due 
to the spatial fluctuations of the local transition temperatures one can find a macroscopic number 
of "ferromagnetic islands" , well separated in space, that spend most of the time in state with a non- 
zero local magnetization which can be either positive or negative. As long as these islands are rare 
(i.e., away from the supposed ferromagnetic transition temperature), they lead to the existence of 
an exponential number of local minima. Moreover, the presence of rare exponentially large islands 
results in the existence of non-analytic (Griffith- like) contributions to the thermodynamic functions 

0- 

An indirect indication that the phase behavior of such systems could be more complicated than 
that described by the renormalization group has been obtained in the framework of the so-called 
non-perturbative renormalization group (RG) approach p. In this latter, the existence of many 
different local minima of the disordered Hamiltonian is taken into account in the form of a replica 
symmetry breaking scheme, and it was eventually found that the renormalization flow leads to the 
strong coupling regime at the finite spatial scale, and not to the expected fixed points. This may 
indicate that something is basically wrong with the supposed (trivial) minimum of the renormalized 
Hamiltonian. 

It has been suggested that when lowering the temperature, the localized ferromagnetic islands 
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become close and strongly interacting, which leads to a transition to the global ferromagnetic state 
0. The solution of the saddle-point equations within the Gaussian variational approximation and 
the replica framework described in the next sections shows that this is not the only possibility. 
We indeed find that upon lowering the temperature the global state of the system can become a 
spin glass before the ferromagnetic state sets in. In this spin-glass state the total magnetization 
remains zero, and there is an effective freezing of local (random) spin configurations (which leads 
to a non-zero value of the spin-glass Edwards- Anderson order parameter). Besides, one finds that 
in this spin-glass state the two-point spin-spin correlation function is described by a temperature 
independent finite correlation length (it is interesting to note that this length coincides with that at 
which the strong coupling regime of the RG approach sets in), while the (spin-glass type) four- 
spin correlation function becomes critical at the spin-glass phase transition point. Finally, when 
further lowering the temperature the global ferromagnetic state eventually sets in via a first-order 
phase transition. 

The existence of an intermediate spin-glass phase in a system where a priori no frustrations, no 
competition of interactions occur is puzzling. This was stressed by Sherrington in a response to 
the perturbative analysis (not within the framework of the replica method) of Ma and Rudnick 
that predicted such a spin-glass phase. Besides the potential flaws associated with the perturbative 
treatments, the problem lies in the fact that it is hard to imagine a disordered ferromagnrt in a 
state where the 4-spin spin-glass susceptibility is larger than the square of the 2-spin ferromagnetic 
susceptibility p| nor in a state with a zero total magnetization and a non-zero spin-glass order 
parameter. This will be discussed later, but, at the level of "hand- waving arguments", one can 
propose a possible interpretation for the presence of a spin-glass phase. The point is that as one 
lowers the temperature (from the paramagnetic phase side) and the ferromagnetic islands become 
close and strongly interacting, there need not be the appearance of a unique infinite (percolating) 
ferromagnetic island. The existence of the spin-glass solution of the saddle-point equations in the 
considered disordered Ising ferromagnet requires that a large number of effectively independent 
spanning ferromagnetic clusters appear in the system. Just below the transition each of the clusters 
is characterized by non-zero value of its own global magnetization (so that within the cluster the 
spins are effectively "frozen"), but the sign of these magnetizations remains random from cluster to 
cluster. This situation manifests itself as the spin-glass state with "frozen" spins and no (averaged 
over clusters) global magnetization. Finally, when the temperature is further decreased, the effective 
interactions among these spanning clusters become strong enough for the system to eventually make 
a "jump" (via a first-order transition) into the ferromagnetic state. In other words, the ferromagnetic 
phase sets in due to a collective locking of the orientations of the clusters magnetisations in the 
same direction. It is easy to understand that this transition must be first-order. Indeed, since at the 
point of the spin-glass to ferromagnetic transition the absolute value of the (randomly directed) 
magnetizations of the spanning ferromagnetic clusters in the spin-glass phase is already finite, 
the value of the global ferromagnetic order parameter resulting from the locking of the various 
orientations in the same direction is itself finite. 

In the next section we present the general formalism in terms of the standard replica approach 
and of the Gaussian variational approximation as applied to the random temperature model. 



In Section 3 we derive all the solutions of the corresponding saddle-point equations, solutions that 
describe the different "ground states" that can exist in the model. It is shown that the spin-glass 
solution discussed above can exist only in dimensions D < 4. Since the conclusions of the present 
study are, to a large extent, only of qualitative nature, we focus on the system in dimension D = 4 
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(the generahzation of the results for dimensions D = {A — e) is given in the Appendix C). We 
obtain the solutions for the paramagnetic, (replica-symmetric) spin-glass and ferromagnetic states, 
and we derive the temperature regions over which these phases are stable as well as the nature 
of the phase transitions separating these these phases. In section 4 the singularity in the spin- 
glass-type four-spin correlation function and in the corresponding susceptibility at the spin-glass 
phase transition is derived. (In Appendix A we give the formal proof that, in the framework of 
the present formalism, no replica symmetry breaking solutions, either continuous or step-like, can 
exist in the spin-glass state.) Finally, in section 5 we discuss the validity of the Gaussian variational 
method; we stress, in particular, that the present approach can only be reliable for finite values of 
the parameter describing the disorder strength. We also suggest a possible scenario for the existence 
of an intermadiate spin-glass phase. 



2 General formalism 

In this paper we study the disordered (random temperature) D-dimensional Ising ferromagnet 
which can be described in the continuous by the following Ginsburg-Landau Hamiltonian: 



H[(P{x);St{x)] 



d^x 



\ {Vmf + \{r- 6r{x)) cl>\x) + \gcl>\x) 



[2.1] 



Here, r = (T — Tf^jTf. ^ 1 is the reduced temperature, and the quenched disorder is described by 
random spatial fluctuations of the local transition temperature 5r(x) whose probability distribution 
is taken to be symmetric and Gaussian: 



P[5r] =poexp(^-i-|rf^a;(5r(a;))2^ 



:2.2) 



where u is the parameter which describes the strength of the disorder and po is an irrelevant 
normalization constant. 

In terms of the standard replica method, the averaged (over quenched disorder) free energy is 
calculated from an the annealed average involving n copies of the same system: 



-(InZ) 



lim —In 

n^O n 



:2.3) 



where (...) denotes the averaging over the random function (5r(x) with the probability distribution 
O), and 



Z-^ = / VbT{x)P\bT\ 



V(f){x) exp {—H [(t>{x); t{x)]) 
is the replica partition function. Simple Gaussian integration over 6t{x) in eq.( p.4|) yields 



(2.4) 



a=l 



V4>a{x) 



exp 



(2.5) 



where 
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is the replica Hamiltonian and 



in i n in 

a=l a=l a, 6=1 



[2.6) 



gab = gSab - u 



[2.7) 



To take into account the possibihty of ferromagnetic ordering in the system we exphcitly intro- 
duce the ferromagnetic order parameter m = (0) by redefining the fields as follows: 



(f)a{x) =m + (fa{x) , (2.8) 
where the new fields ipa{x) describe the spatial fluctuations with zero mean. By substituting (p.8|) 



into eq.(2^) for the replica Hamiltonian, one finds 



i/(")[(^a(x);m] = Vn{Km^ + ^gm^) + J d^x 



2 n 1 " 1 " 

1^ J2 (VV'a)^ + 9^ E V'a + T J2 9ab^l^l + 
^ a=l a=\ ^ a,b=l 



+ ^9'^'^ V^a + J2 gabVaVb + Tm ^ + m ^ gab^a^l + nt' g 

a=l a, 6=1 a=l a, 6=1 a=l 



(2.9) 



where V is the volume of the system. Note that the limit — > 0, that must formally be taken in 
the final results, allows us to omit all terms of order rt^ in the above expression (and in further 
calculations). 

The idea of the Gaussian variational approach is to approximate the fluctuations of the fields 
iPa{x) in the above eg. (|2.9| ) by the Gaussian trial Hamiltonian 



V 



n 



Gjip)iPaiPWbi-p) 



(2.10) 



2 Jm<i (27r)^ 

where the correlation functions Gabip) = {Va{p)'^b{—p)) are considered as variational parameters 
The replica partition function can be represented as follows: 



n 

a=l 



Vipa{x) 



exp{-i7(")[y,,(x)] - (//(") [y.. (x) ; m] - H^-\ipa{x)])] 



(2.11) 



and in the first-order cumulant approximation in the difference (if*^") — if^"-*) one finds 



~ exp 



= exp(-nK/[m;G]) (2.12) 



where {{...))g denotes the averaging with the Gaussian weight, eq.( p.lO|) ; / [m; G] is the density of 
free energy that depends on the order parameter m and on the trial correlation functions Gab{p)'- 
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/ [m; G] = ^ / ^Tr HG-\p)) + - H^^^)), . (2.13) 

J\p\<i [Ztt)^ nV ^ 

Since the above free energy density is an upper bound of the exact rephca free energy density, 
the variational parameters m and Gab{p) can be determined by minimization of eg. ( P.13|) . One 
should however keep in mind the oddities related to the limit n — 0, in particular the fact that the 
number of parameters can turn negative for n < 1 (see section 3.3)^. Inserting eqs.( p.9|) and (|2.1CI|) 
into eq.( p.l3| ) leads to 



/ [m; G] = -— / T^Tr ln(G(j9)) + -rm^ + -gm^ + 
^ ' ^ 2n J\p\<i (27r)^ ^ 2 4^ 

1 r d^p " 1 " 



a=l a, 6=1 

n \ ^ 

+ 7r9^^Y.('^l(^))9 + J2 9ab{Va{x)<^b{x))g . (2.14) 

a=l ^ a,b=l 

Above and in what follows we omit irrelevant constant terms. For the Gaussian averages of the 
fluctuating fields one has 



H<i (2^0^ 
d^p 



Gab{p) = [Gab] 
Gaa{p) = [Gaa 



I\P\<1 (2vr)^ 

- {ipl{x))Mi^))9 + 2{^a{x)Mx))l ^ [Gaa] [Gbb] + 2 [Gab] 

where we have introduced the notation 



{'^lix)'^bix))c 



I 

J\v\ 



d'^p 



A{p) ^ [A] 



(2.15) 

(2.16) 
(2.17) 



(2.18) 



hp\<i (Stt)^' 

for an arbitrary function A{p). Taking into account that the diagonal elements of the matrix G 
must be independent of the replica index, Gaa = G we find the following expression for the free 
energy density: 



f[m-G] = -i-Tr [ln(G)] + ^rm2 + l^m^ + ![(/ + 



^ As noted by Mezard and Parisi in their replica field theory for random manifolds ||T^ , the Gaussian variational 
method becomes exact when the number N of components of the fields ipa goes to infinity (here, N = 1). To apply 
this remark to the present case, one must generalize the non-Gaussian term appearing in the replica Hamiltonian, 
eq.(2.6), to X]ab=i(5<ifc/A^)[0a^fc + 2(0^ • 0;,)^]. Notice that with this latter term the Hamiltonian does not 
correspond to the replica-space formulation of the random temperature 0{N) model. 
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G 



2 1 



2 1 



+ -{g-u)[G\ + ^ E 9ah [GahY + 



1 r~-| r~-| 1 " 

+ -gm^ G + m^{g - u) G + -m^ V Qab [Gab] 

2 L J L J 



(2.19) 



The correlation functions Gabip) and the order parameter m are then determined by the follow- 
ing saddle-point equations: 



8f 



6G{p) 

Sf 
SGM 

51 
5m 



{a^h) 



By using the explicit expression of the free energy density, eq.( p.l9| ), one obtains 



(2.20) 
(2.21) 
(2.22) 



Gabip) = (P^ + 'T)5ab + 9 G 5ab + ^Qab [Gab] + Q^^^ab + '2^^^ Qab 



(2.23) 



m 



r + gm^ + {:ig-2u) G + -Y,gab[G, 



2 A 



^ab\ 



a-^b 



(2.24) 



According to eq.( p.23| ) one finds that the trial correlation function has the following structure: 



Gab (P) = (P + ^)5ab + f^ab 



(2.25) 



where the matrix fiab is defined by 



l^ab = [9 



G 



+ giv?) Sab + 2gab [Gab] + 2gabm^ ■ 



(2.26) 



For finding explicit solutions of this equation one needs to make an assumption about the replica 
structure of the matrix fiab- In what follows we assume that this matrix is replica symmetric; 
in Appendix A we give the formal proof that eq. (|2.26| ) has no solutions with the Parisi replica 
symmetry breaking structure for the matrix fiab- The replica symmetric ansatz implies that the 
matrix fiab is defined by only two parameters, 



/iab = (/i + /x)5afe = I (2.27) 

For the corresponding replica symmetric correlation function, defined by eq. (|2.25 ), we find 
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Sab + 



= G'c(p;A) 5afe + /i(G'c(p;A))^ 
where the so-called "connected" part of the correlation function is given by 



Gcip; A) 



1 



(2.28) 



(2.29) 



p2 + A 

and instead of fl (defined in eq.( f2.27D ) we have introduced a physically motivated "mass" parameter 
A = T + jl + that determines the value of the correlation length {Rc ~ X^^^"^ in the present 
approximation) . 



Note that according to eg. ( |2.28| ) the parameter n is related to the value of the spin-glass 
Edwards-Anderson (EA) order parameter, since 



n— >0 

Thus, to be physically meaningful fi must be non-negative. 



(2.30) 



By using eqs.( |2.27| ), ( |2.28| ) and (p.7|), the corresponding saddle-point equations for the parame- 
ters m, A and fj, can be obtained from eqs. (|2.26|) and ( |2.24|) as 



\ = T + fi+{3g-2u) ([Gc] + /i Gl) + {3g- 2u)m 



fl = 2ufi GI + 2um 



m{T + gm^ + {3g - 2u) (^[G^] + n ) + 2ufi 
The resulting the free energy density is then given by 



G 



(2.31) 
(2.32) 
(2.33) 



/(m,A,/.) = -l[ln(G,)] + ^(r-A)([GJ+/i[G,2])+^(3^7-2n)([a]+/x 



2 

1 2 
+ -Ufl 



G„ 



+ 



G 



+ -{3g-2u)m'[[G,]+ii 



G„ 



+ uixm 



G 



+ 



1 2 1 4 

H — Tin + -qm 

2 4 



(2.34) 



3 Phase diagram in D = 4 



In this Section we study all possible solutions of the saddle-point equations ( p.31|) -( p733D . As usual, it 



is assumed that the non-Gaussian coupling parameters g and u of the original replica Hamiltonian, 
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eq.( p.6|) , are small: g <ti 1, u 'ti 1. Besides, we consider u g, which qualitatively corresponds to the 
situation of "finite disorder strength" since then the parameter u describing the disorder strength of 
is of the same order as the coupling parameter g of the pure system. As will be discussed in section 
5, the analysis of the validity of the present (first-order) Gaussian approximation shows that it may 
give reasonable results provided the ratio u/g stays within certain numerical bounds (see Section 
5). For the moment, however, it is sufficient to assume that u < ^g (see below). Finally, since we 
are only interested in the large-scale (continuous limit) properties of the system we consider the 
region of parameter space where the mass A of the connected correlation function is also small: 
A < 1. 

To simplify the algebra and for a qualitative presentation of the phase diagram it is convenient 
to consider first the solutions of the saddle-point equations in dimension D = 4. Generalization of 
the results for dimensions below four will be given in Section 5 (it will also be shown that for D > 4 
the spin-glass solution does not exist). 

For D = 4 and for A ^ 1 one has 



Gl 



|p|<i (27r)>2 + A V ^A- 



b|<i (27r)4 (p2 + A)^ 



C In 



(3.1) 
(3.2) 



where C = l/167r^. 



3.1 Paramagnetic solution 

In the paramagnetic state, the ferromagnetic and spin-glass order parameters are both zero (m 
= 0) and there is only one saddle-point equation (|2.31|) for the mass parameter A 



A = r + (3^ - 2u) [Gc] . (3.3) 
Using eq.( p.l|) leads to the following equation: 

A + C{3g - 2n)Aln = r + G{3g - 2u) (3.4) 

which provides the dependence A = A(r). The solution of this equation makes physical sense 
only for A > 0, and therefore this condition defines the limit of existence of the paramagnetic 
phase. Provided 3g > 2u the above equation yields positive (physical) solutions for A(r) only for 
temperatures such that 



T>Tc = -G{3g - 2u) . (3.5) 

If one finds that the ferromagnetic solution appears just below Tc, then the temperature r = Tc 
would correspond to the paramagnetic- ferromagnetic phase transition point. However, it will be 
shown below that this is not the case. In fact, below a certain temperature Tgg > Tc a spin-glass 



9 



solution (with /i 7^ 0) appears, and at temperatures r < Tgg it is the spin-glass state that turns out 
to be stable, while the paramagnetic one becomes unstable. 



3.2 Spin-glass solution 

The spin-glass state is defined by two saddle-point equations (p.31|) and ( p. 321) : 



T + fi+{3g-2u)([Gc]+fi 



Gl 



(3.6) 



2u/i 



(3.7) 



which define two non-zero (positive) order parameters: ^{t) and A(r). From the last equation one 
immediately finds that for 7^ the mass parameter A becomes temperature independent, A = Aq, 
and the value of Aq is defined by the condition 



1 

2u 



Correspondingly, eq. (p.6| ) yields the following solution for the spin-glass order parameter: 

2?/ 

/i(^) = TT (^o - {^9 - 2u) [GJ - r) . 
3^ 

By making use of eqs.(^.lD and ( |3.2| ) one can find the solutions for Ao and /i(r) explicitly: 



(3. 



(3.9) 



where 



Ao 



exp 



2Gu 



- 1 



T 



(3.10) 
(3.11) 



Tsg = Ao - {"ig - 2u) [Gc] =Tc + 



^ + G(3, 
2u ^ ^ 



2u) Ao > Tc 



(3.12) 



and T(. is the putative paramagnetic critical point discussed above. The solution for fi > appears 
(i.e. becomes physical) only for r < Tgg, and therefore the point Tgg can be associated with the 
spin-glass phase transition temperature. Note that according to eqs.( |3.1^ ) and ( |3.3| ), the value of 
A in the paramagnetic phase at r = Tsg is equal to Ao (for r > Tsg, A(r) > Ao and for r < Tsg, 
A(r) < Ao). Since nijsg) = whether one comes from the paramagnetic or the spin-glass phase, the 
transition into the spin-glass phase is clearly continuous. 

In Appendix B we present a detailed study of the stability of the spin-glass and the paramagnetic 
solutions obtained above. It is shown there that for r > Tgg the only stable state of the system is 
paramagnetic, while for r < Tsg the paramagnetic solution becomes unstable and the stable state 
of the system is the spin-glass phase. 
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3.3 Ferromagnetic solution 



We finally consider the ferromagnetic solution of the saddle-point equations ( |2.31| )-( p.33| ) in which 
all three parameters A, fi and m are non-zero. After some simple algebra we find 



m 



u 



X 



/i 



gl-2u 

where the parameter A(r) is obtained from the following equation: 



(3.13) 
(3.14) 



A = r + (3^ - 2u) [G] 



A 



2 1 - 2m [G2] 

Substituting eqs.( |3.1| ) and ( |3.2| ) into the above equation gives 



(3.15) 



A 



| + C(39-2.)-C(39-2.)lnA__i_ 1 



T + C{3g-2u) 



(3.16) 



A simple analysis shows that upon lowering the temperature r a solution of this equation appears 
for the first time below a temperature r*. This solution has a finite (non-zero) value A,, at r = r^,, 
which indicates that the phase transition into the ferromagnetic state is first order. To the leading 
order in -C 1 and in m -C 1 (and for ^f/w ~ 1) one finds 



r, ~ -C{3g - 2u) 



exp 



1 



2Cu 



Tr — 



3e 



Ao < Tc 



(3.17) 



and 



A(r = r*) = A* 



exp 



2Cu 



eXo ■ 



(3.18) 



By inserting the above value into eqs.( 3.13 ) and ( |3.14 ), we find the corresponding values of the 
ferromagnetic and the spin-glass order parameters: 



Straightforward calculations similar to those of Section 3.3 show that the ferromagnetic solution 
defined by eqs. (|3.13|) - (|3.15|) is stable at all temperatures r < r*. Thus, below both the spin-glass 



and the ferromagnetic solutions are (locally) stable (this is the standard situation for first-order 
phase transitions). To determine which of these two states is the global minimum of the free energy 
at a given temperature we have to compare the corresponding values of their free energies. 
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3.4 First-order phase transition between spin-glass and ferromagnetic 
states 



Substituting the spin-glass solution, eqs.(p.8|)-( |3.12| ), into eq.( p.34| ) provides the value of the free 
energy density of the spin-glass state (in the leading order in g^u <^ 1): 



T -Tr 



where 



(3.21) 



u 



u 



1 



/o(r) = -C{?>g - 2u) + -r - -r' . 
6g 3g 12 

On the other hand, for the ferromagnetic solution, eqs.(|3.13|)-(|3.16|), one gets 



(3.22) 



//( 



/o(r) + ^(l + 2CMlnA) (r - r. 



C'^u{3g - 2u) 



Qg 



+C 



Qg 



4 



1 
24 g 



1 



where the value of A(r) is given by eq.(|3.16|). Let us redefine 



(3.23) 



A(. 



Xox{t) 



(3.24) 



-Xot (3.25) 
By making the above change of variables in the saddle-point equation ( p.l6|) , one obtains 



X 



| + C(3«-2„)(l-.u.)-^i- 



t . 



(3.26) 



Therefore, to the leading order in u,g <^ 1 the value of the parameter x as the function of the 
reduced temperature t is defined by the following equation: 

4^T^ ^ ^ ■ (3-27) 

Assuming that at the point of the phase transition, i.e., when ff = fgg, the value of the parameter 
x{t) is of order one, we find for the difference of the free energies, eqs.( p.21"D ,( |3.23|) (in the leading 
order in u, g <^ 1 and Ao ^ 1), 



Cu 1 

ff — fgg ~ — A^t (1 — X — x\nx) H A^x^ . 

Ji Jsg o K ' Sg ° 

Thus, the transition point (// — fgg = 0) is defined by the following equation: 



(3.28) 



12 



c 



ut{l~x — xlnx) 



-X 



(3.29) 



Combining eqs.( |3.27|) and ( |3.29| ), we finally derive the equation for the parameter x at the phase 
transition point, 



3 

-X 

2 



xlnx = 1 



(3.30) 



This equation has a unique solution x = x/~l(x/>l). Substituting x/ into eqs.( p.27D and (|3.25|) 
one obtains the temperature of the (first-order) phase transition between the spin-glass and the 
ferromagnetic phases: 



Tr — 



Xf 



4:Cu\nxf ° ' 



(3.31) 



which is less than r*. 



4 Singularities at the spin-glass phase transition 

Within the Gaussian variational approximation the (connected) correlation functions can be ob- 
tained by adding source terms to the replica Hamiltonian, eq. (|2.6|), and by approximating the 



free energy functional to the first-order cumulant as in eq. (|2.12|) . This latter then generates the 
(connected) correlation functions that are obtained by functional differentiation with respect to the 
source terms. By introducing source terms linearly coupled to the fields (j)a{x), one derives the usual 
correlation functions, whose expression (for the two-point functions) coincides with that given in 
section 2. To study the singularity at the spin-glass transition it is more convenient to introduce 
source terms linearly coupled to the composite operators |0a(x)0;,(x). This leads us to the following 
replica Hamiltonian: 

[0,(x); A(x)] = [0,(x); 0] - - / t/^x ^ A„,(x)0,(x)0,(x) (4.1) 

a,b=l 

where H^'^^ [0a(x);O] is given by eq. ( p.6| ) and, because we are ultimately interested in the replica- 
symmetric solution, the source term is taken as 

A,,(x) = (A(x) + A(x)) 5ab - A(x) . (4.2) 

Since we only study here the paramagnetic and the spin-glass phases, we can set m = 0. A 
Gaussian trial Hamiltonian is chosen as before, but the presence of space-dependent source terms 
breaks the translational invariance and requires trial Green functions that depend on two space 
points: 



ab=l 
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where 



y dxdx 



-; + r\ 6ab + ^ab{x) 



5{x — x') 



(4.4) 



which simply generahzes eq. (|2.25| ). For the rephca-symmetric solution, the variational parameter 
fiab{x) can be written as 



f^abix) = {n{x) + Sab - K^) , (4.5) 

whereas the Green functions can be written as 

G:t' = GfSab + Gf (4.6) 

where Gc and Gd represent, as usual in the presence of quenched disorder, the "connected" and 
"disconnected" parts and can be expressed in terms of r, fi{x) and fi{x) by inverting eq.( f4.4|) . 

One can then follow the procedure used in section 2: the first-order cumulant approximation in 
the deviation (H^"-^ —H^^^) provides an upper bound for the free energy functional, and minimizing 
with respect to the trial Green function elements leads to saddle-point equations that in the limit 
n — >^ for a replica-symmetric scheme reduce to 



li{x) + A(x) 



(3(7 - 2u) {GT + Gf 



(4.7) 



/i(x) + A(x) 



2uGf 



When considered at the saddle-point characterized by the above equations, the variational free 



energy F 



Afx), Afx) 



can be used as generating functional for the correlation functions of the 



composite operators |0a(a;)0b(x). More precisely, one has 



5F 
5K{x) 

6F 
6A{x) 



-o lim(0a(a;)</>a(a;)) = --(G; 



+ GT) 



- \im{(f)a{x)(j)b{x)) 

2 n—*0 



[aj^b) 



> d 



(4.9) 
(4.10) 



In the limit A = A = the above equations reduce to the expressions already given in section 2, in 
particular eq.( [4.10D reduces to eq. ( [^.3U| ) for the spin-glass order parameter: indeed, in the absence 
of source terms translational invariance is recovered and Gf = [Gdip)] = [G^c(p)]- 

The second functional derivatives of F provide the wanted two-points, four-field correlation 
functions, namely: 



6^F 



6A{x)6A{x') 



1 1 " 

a,b=l 



(4.11) 
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r2ri 11" /" \ 

sr2 rp 11"/" \/" \ 

where ((...))c denotes a cumulant average for the composite operators. 

We are interested in the case A = A = 0, case in which translational invariance is recovered and 
the two-point functions are diagonal in momentum space. We then define the 2x2 momentum- 
dependent susceptibihty matrix xiv) by 



Xii(p) 



ST 



6A{-p)6A{p) 
6^F 



6^F 



6A{-p)6A{p) 



A=A=o ' ^''^^^ ~SAi-p)5Aip) 
A=A=o ' ^""^^^ ^ SAi-p)6A{p) 



A=A=0 



A=A=0 



(4.14) 



The expression of these p-dependent susceptibihties in terms of the correlation functions follow from 
eqs. (|4.11|) - (|4.13| ). By combining the above definitions, the expression of the free energy functional at 
the first-order cumulant approximation and the saddle-point equations, eq. ( |4.71) -( |4.8| ), one obtains 
after some algebra the following expressions for the p-dependent susceptibility matrix: 



xip) 



( 



\ 



1 



(3s-2«) 









2u 



yr^{p)-\ 



(4.15) 



where I is the identity matrix and the p-dependent matrix M(p) is defined by 



M(p) 



/ l + (3^-2M)(/(p)+2/iJ(p)) 2(3^ - 2M)/iJ(p) 



4m/xJ(p) 



1-2m(/(p)-2Mp)) 



(4.16) 



with 



^^^^ /p|<i (27r)«(k2 + A)((k + p)' + A) ^^'^^^ 

Up) = JJ_M = f ^ I (4 18) 

2 d\ 4|<i(27r)^(k2 + A)2((k + p)VA) ' ^ ^ 



and A and fi are given by eqs. (|3.6|) and ( |3.7| ), respectively. From the above equations one immediately 
obtains the expression for the susceptibility matrix in the paramagnetic phase (where /i = 0): 
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I Hp) 

I l+{3g-2«)/(p) 



V ^(p) / 

\ ^ l-2«/(p) / 



(4.19) 

from which one derives that the susceptibihty X22(p = 0) = {d'^F/dA'^) Ii^^^^^q-j diverges when 



1 — 2m/(0) = 0, i.e., by using eq. (|4.17|) , when [G^] = l/2u. This point is precisely attained at the 
transition to the spin-glass state, t = Tgg, where, for D = 4, Tsg is given by eq.( p.l2| ) (A is then equal 
to Ao given by eq. (|3.1CI|) ). To derive the critical behavior of X22{p) when approaching the spin-glass 



transition from above, we use the small-p expansion of I{p) in D = 4, 

1\ .\ C2 2 



lip) ^ C ( In ( - 1 - 1 1 - (4.20) 



where C2 > 0. After defining A = Ao + 6\, 6\ 0+, this gives 



C2 

Using eqs.(|3^), ( p.lO| ) and ( |3.12|) one finds that when r - 



sgi 



3<7 

which finally leads to 



A. 



27/ 

SX~^{t-t.,) , (4.22) 



Consider now the spin-glass phase (in D = A). One then has 2m/(0) = 1, A = Aq, and /i = 
1^(^89 ~ t), so that for small p, 

I{p) ^ (4.24) 
lu Ao 

J(P) ^ ^ - • (4.25) 
The determinant of the matrix M(p), eq. ( |4.16| ), can now be expressed at the leading order as 

detmv)) - ^ + |^(-^. - -)) (4-26) 
so that when r — > r^^, the susceptibilities Xiiip) stay Xi2(p) are finite, whereas 

16 



One concludes from the above formulas that the susceptibility associated with the external field 
A that couples to the Edwards- Anderson order parameter g, i.e., X22(p) = sA{-p)5A{p) ' diverges 
when the critical point Tgg is approached both from above and below as |r — Tsg\^^ (in D = A), 
whereas the associated correlation length (that characterizes the long-distance behavior of the two- 
point composite- field correlation function, eq.( [4.13|) ) diverges as |r — T^gl"^/^ (in D = A). The 
corresponding critical exponents, 7 = 1, = 1/2, are thus classical. 



5 Discussion 



5.1 Dimensions other than D = 4 

The situation in dimensions D > A is quite simple. There, the value of the integral 

^ r d^p 1 

~J\p\<i {2n)Dp2 + X ^ ■ > 

remains finite (not diverging) in the limit A ^ 0. Therefore when m <^ 1 the only solution of the 
saddle-point equation ( p.7|) for the spin-glass order parameter is trivial, yU = 0, so that there is no 
spin-glass solution in dimensions D > A. Thus, one recovers in this case the standard scenario: upon 
lowering the temperature, the only phase transition that takes place in the system is a second-order 
phase transition from the paramagnetic to the ferromagnetic phase, and this phase transition is 
described by the Gaussian theory. 




On the other hand, the phase diagram in dimensions D < A turns out to be similar to that in 
D = 4 at a qualitative level. The integral, eq.( |5.1|) , diverges in the limit A ^ 0, and, therefore, 
when M ^ 1 there is always a spin-glass solution 7^ of the saddle-point equation (|3.?| ). Thus, 
in this case one recovers within the Gaussian variational approximation a phase diagram similar to 
that in dimension Z) = 4, where the paramagnetic phase is separated from the ferromagnetic one 
by an intermediate spin-glass phase. 



5.2 VaUdity of the Gaussian variational approximation 

Since the main result of the present study, namely, the existence of a spin-glass phase separating the 
paramagnetic and ferromagnetic phases, is in apparent contradiction with the generally accepted 
view on the phase diagram of the disordered Ising ferromagnet, the limits of validity of the Gaussian 
variational approximation used in this paper requires a detailed study. 

It is well known that the Gaussian variational approach becomes exact when the number of spin 
components tends to infinity (see footnote on page 5). Otherwise (in particular, the Ising model 
is very far from this limit) it is not more than an approximation characterized by certain bounds 
of validity (if any), and ignoring these bounds may just lead to wrong conclusions. The typical 
example is the pure Ising model in dimensions D < A: there, one can easily check (using eqs. (|2.31|) 
and (|2.33|) with m = /i = 0) that according to the Gaussian variational approximation the transition 
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from the paramagnetic to the ferromagnetic phase turns out to be first-order, which is of course 
incorrect. 

The vahdity of the first-order cumulant approximation in the deviation (iJ*^"-* — i/^"^), described 
in section 2, can be checked by estimating the contribution from the higher order terms. In the 
vicinity of the critical temperature Tc of the (supposed) paramagnetic-ferromagnetic phase transi- 
tion, for \T — Tc\/Tc = r ^ 1, a systematic account of these contributions can be done in terms of 
the usual renormalization group (RG) procedure, which yields an effective scale dependence of the 
renormalized non-Gaussian interaction parameters g and u. In dimensions D = 4, in the so-called 
one-loop approximation, the scale evolution of g and u is described by the following well-known 



RG equations (see e.g. m, M, [12 



-g{0 = -QC {3g -Au)g + (g^- g\ gu^; u^) (5.2) 

-^uiO = -4C (Sg -Au)u + {g'; g\ gu'; u') (5.3) 

where, as usual, C = l/lGvr^, and ^ = InL is the standard RG rescaling parameter which is equal 
to the logarithm of the spatial scale. According to these equations one can conclude that the 
renormalization of the parameters g and u remains irrelevant (so that the higher order terms of the 
perturbation theory are not important) and that the theory remains Gaussian at scales bounded 
by the conditions 



6C\3go~4uo\go^<^go (5.4) 

AC\3go - Auo\uo^ <^;: uo (5.5) 

where g^ = g{^ = 0) and mq = = 0) are the "bare" (microscopic) values of the parameters g 
and u. These two conditions are satisfied when 



6C\3go-Au 







or, in terms of the spatial scale L, until 



L<^L,r^ exp{ ^ , , } • (5.7) 

Since the temperature and the spatial scales are related by |r| ~ L^^, one finds that the higher-order 
non-Gaussian corrections remain irrelevant only at temperatures not too close to the (supposed) 
paramagnetic-ferromagnetic critical point: 

|r| > ~ exp{- ^ , , } • (5-8) 

3C\3go - Auo\ 
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Note that in the special case uq = ^go, the higher-order terms of the RG eqs.(p.2|)-( ^.3|) come 
into play, and instead of the condition (|5.8|) one eventually obtains: ~ exp{ — (const) / g^} . Note 
also that in the case of the pure system, uq = 0, eq.( |5.8|) yields the well known temperature scale 
~ exp{ — 167r^/9(7o}, such that when |r| ^ the behavior of the system is effectively Gaussian, 
while in the close vicinity of the critical point, for |r| <^ r*, the non-Gaussian fluctuations become 
dominant; in particular, this shows that the first-order Gaussian cumulant approximation breaks 
down for |r| ^ r^, (it is at this temperature scale that this approximation wrongly predicts the 
first-order nature of the phase transition). 

According to the calculations of section 3, the spin-glass phase separating the paramagnetic 
and ferromagnetic phases exists in a temperature interval of order of Ao — exp{ — j^} around 
(where is the critical temperature of the putative paramagnetic- ferromagnetic phase transition). 
Since the approximation used in the present approach is valid only outside the temperature given 
by eq.(|5.8|), we need this "dangerous" temperature interval to be well inside the spin-glass phase. 



If the above condition is satisfied, the spin-glass state (obtained in terms of the present Gaus- 
sian variational approximation) would set in when lowering the temperature well before the non- 
Gaussian critical fluctuations of the fields 4>a{x) (given by the higher-order terms of the perturbation 
theory) become relevant. 

It should be stressed, however, that this does not guarantee that the critical behavior at the 
spin-glass phase transition itself is correctly described by the Gaussian theory (as it is derived in 
section 4), since for an adequate description of this phase transition one needs to take into account 
critical fluctuations of the spin-glass order parameter qab{x) = (j)a{x)(f)b{x) and not just those of 
the fields (pa{x). Actually, and although less likely, one can not exclude that the fluctuations of the 
composite field qab{x) will even destroy the spin-glass phase itself. This question is left for future 
investigations. 

In addition, when further lowering the temperature, the transition point from the spin-glass 
to the ferromagnetic one is also separated from by an interval of the order of Aq. Therefore, 
provided the condition (|5.9|) is satisfied, one can conclude that the first-order nature of this phase 
transition corresponds indeed to the proper physical phenomenon, and is not just an artifact of the 
method used. 

Since all the calculations of the present study are done for (? <^ 1 and m <^ 1, one can easily see 
that the condition ( |5.9| ) is satisfied provided 

10 Oo 14 , , 

so that the parameter u describing the disorder strength must be of the same order as the interaction 
parameter g of the pure system. It is in this sense that we characterize such system as having a 
finite strength of disorder (although both u and g are kept small). 

Thus, we have found at least a limited region of the parameters, defined by eq. ( ^.10|) , where the 



approximation used in this paper appears to be reasonable and for which an intermediate spin-glass 



19 



phase could exist. Presumably, it could exist also in a wider region (in particular for large values of u 
compared to g), where, however, it cannot be studied by the present method. Unfortunately, in the 
framework of the present investigation we cannot answer the question of whether the (exponentially 
narrow) intermediate spin-glass phase continues to exist in the limit -u ^ or whether there exists 
a critical value Uc such that the spin-glass phase appears only for u > Uc- 



At a qualitative level, the situation in dimensions D = 4 — e(e^l) turns out to be similar 
to that m D = A (see Appendix C): we find again the same restriction on the parameters g and u 
as in eq. (|5.10|) , but in addition we obtain that the value of the parameter e = A — D must remain 
small. This shows, on one hand, that the results obtained in the present study are not unique to 
the dimension D = A and can be continued down to lower dimensions; on the other hand, due to 
the restriction e <C 1, we cannot guarantee that they would survive down to the physical dimension 
D = 2, (where, for sure, the present approximation does not work). 



5.3 Conclusion 

The existence of a spin-glass phase in a "dilute ferromagnet" is a priori a puzzle In the 

absence of competing antiferromagnetic interactions and of the associated frustration, how can the 
system ever be found in a state with frozen magnetizations but without global magnetization?^ 
Indeed, for a spin-glass phase as predicted by the Gaussian variational approximation to exist in a 
disordered Ising models, one must have 

(i) non-zero "frozen" local magnetizations in an extensive part of the volume, so that the 
Edwards-Anderson order parameter is non-zero; 

(ii) a total magnetization equals to zero. 

Truly frozen local magnetizations require symmetry breaking and can exist only on a percolating 
cluster that diverges in the thermodynamic limit. If L is the linear size of the whole system (ulti- 
mately, L oo), the size of clusters with non-zero magnetization must scale as L'^ with D > d > 1 
(for Ising spins). On the other hand, on each cluster the magnetization can be either positive or 
negative. One possible scenario to explain the occurrence of the spin-glass phase with zero total 
magnetization is thus as follows. 

We envisage that there exist a global temperature at wchich a large number M of percolating 
clusters of fractal dimension d < D acquire a non-zero magnetization (when L oo), while being 
essentially decoupled from each other. These clusters are formed from the " islands" characterized by 
a predominantly negative value of the local temperature (r — 5t) and are separated from each other 
by regions of high local temperature. In zero external magnetic fieldQ the sign of the magnetization 
on each cluster is then randomly determined (and subsequently frozen when the temperature is 

^ The fact that at and around the spin-glass transition, the 4-spin spin-glass susceptibility has to be larger than 
the square of the 2-spin ferromagnetic susceptibility may not be as severe a problem: the presence of a large number 
of " ferromagnetic islands" and clusters invalidates the inequalities obtained via a perturbative analysis |^ 

If one imposes a small external magnetic field, all cluster magnetizations will have the same sign, i.e., the sign 
of the external field. Then, in the limit of vanishingly small field, the overall spontaneous magnetization does not 
vanish. The situation, however, is different if one considers from the beginning the case without any applied magnetic 
field. 
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decreased). Since M is large, there will be essentially (up to a relative factor of order \/\fM) as 
many positive as negative clusters. So if M scales with L as M ~ L^~'^ (and d < D), the total 
magnetization, averaged over all clusters, is zero whereas the Edwards-Anderson order parameter 
is non-zero: the system is then in a spin-glass phase. It must be stressed that the appearance 
of the large number of percolating clusters can not be envisaged as a purely geometric, disorder- 
controlled phenomenon, thermal fluctuations come into play for helping the islands to coalesce. This 
is necessary in order to have M scale as L^~'^ since in usual percolation the number of incipient 
spanning clusters at the percolation threshold can be larger but in D < Q does not diverge with L as 
a power law [0]lt must be stressed that the appearance of the large number of percolating clusters 
can not be envisaged as a purely geometric, disorder-controlled phenomenon, thermal fluctuations 
come into play for helping the islands to coalesce. This is necessary in order to have M scale as L^~'^ 
since in usual percolation the number of incipient spanning clusters at the percolation threshold 
can be larger but in D < 6 does not diverge with L as a power law [|l^ 

To settle the question of the existence of a spin-glass phase in the random temperature Ising 
ferromagnet, especially in D = 3, further investigation is clearly needed. A potential line of inves- 
tigation is a renormalization group analysis that treats on equal footing the primary fields (paix) 
and the composite fields 0a(x)06(x) that are needed to describe the spin-glass phase. This would 
clarify the problem of the upper critical dimension of the paramagnetic to spin-glass phase tran- 
sition (recall that in our study D = 4 is the critical dimension above which the spin-glass phase 
ceases to exist, it is not necessarily the upper critical dimension). An interesting problem that is 
worth pursuing is the connection to the phase behavior of the random field Ising model. It has 
been suggested fl^, |]15| that this latter system possesses an intermediate spin-glass-like phase, and 
that a potential source for the presence of this phase is the generation in higher orders of perturba- 
tion theory of "attractive" terms of the type u(f)'^(f)l |]15|, i.e., terms similar to those present in the 
random temperature replica Hamiltonian. Although the kind of symmetry breaking leading to the 
spin-glass phase in the random temperature Ising model is not possible in the random field models 
(within the replica formulation, random fields generate a nonzero field A that is linearly coupled 
to the composite operator (j)a{x)4>b{x)), a comparative study may prove fruitful. Finally, one should 
stress that absence of the replica symmetry breaking found within the present approach could turn 
out to be an artifact of the Gaussian variational approximation because the corresponding saddle- 
point equations (see Appendix A) are in fact rather "fragile" with respect to the introduction of 
higher-order terms of the perturbation theory. 
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6 Appendix A 



Consider again the saddle-point equations (|2.26|) in the spin-glass state (m = 0) 



where 



and G{p) = Ga=b{p)- 



l^ab = 9 



G 



Sab + '^9ab [G. 



ab\ 



Gabip) = {P" + T)5ab + ^lab 



(6.1) 



(6.2) 



6.1 Continuous RSB 

In the case of a Parisi-like continuous RSB in the limit n ^ the matrices Gab and ^ab are 
parametrized by their diagonal elements G, jl and by off-diagonal functions G{p; x), fi{p; x) defined 



on the interval x G [0, 1] |TT[: 



Gab{p) 



G{p)] G{p;x] 



(6.3) 



f^abip) (^(p); fi{p;x)) (6.4) 
By using the Parisi algebra for inverting hierarchical matrices [|1^] one derives from eq.( |6.2|) 



where 



G{p; x) 



/i(0) 



(p^ + r + /i - /i)2 Jo (p^ + A(y))2 



dy- 



(6.5) 



Substituting eq. 
function li{x): 



Kv) = r + fi-fi- yfi{y) + / dzfi{z) 

Jo 



/i 



dxjji^x) 



(6.6) 
(6.7) 



^(0) = ^{x = 0) . (6.8) 
into eq.( |6.1| ) (for a ^ b), we find the following equation for the unknown 



fj,{x) = 2u 



d^p 



+ 2u 



d^p fi'{y) 

dy- 



l\p\<i (27r)^ (j92 + r + ^ - 71)2 ' (27r)^ Jo ^ {p^ + \{y)y ' 

Differentiating both sides of this equation with respect to x then leads to 



(6.9) 



22 



2u 



. 



/h<i (27r)^(p2 + A(x))2 
Thus, either one has fi'{x) = 0, i.e., = const (independent of x), or 



(6.10) 



2u 



1 = 



/H<i(27r)^(p2 + A(a;))2 ^ " ^^'^^^ 

which also yields fi{x) = const. In both cases the solution is replica-symmetric (or step-like). This 
proves that the saddle-point equation (|6.1|) cannot have solutions with continuous RSB. 

Nevertheless, the above proof does not exclude the possibility of other types of solutions with 
step-like RSB. 



6.2 One-step RSB 

We consider now a one-step RSB ansatz for the replica matrix Gabip)'- it is defined in terms of one 
diagonal element g{p), two off-diagonal elements gi{p) and go{p), and the coordinate of the step Xo'. 



Gabip) 



g{p) ; for a = 6 

91{P) ; for /(!;)= /(A 



(6.12) 



where I{x) is the integer valued function which is equal to the smallest integer larger than or equal 
to X. Substituting the above equation into the general expression, eq. ( p.lQj) , (with m = 0), one 
finds the following explicit expression for the free energy density: 



F [g{p) \ gi{p) \ gQ{p)\Xo] 



2x 



[In {9-91+ Xo{gi - go))] - — ]\ki{g- gi] 



Xq 



9o 



'^{9-91 + Xoigi - go)) 



+ 2 



(/ + r)g 



+ 



+\i^9 - 2m) [gf + ^m(1 - Xo) [gif + ^uxo [go]'^ 



(6.13) 



where, as usual, the symbol [(...)] denotes the integration over p, see eq. ( |2.18| ). Variation of this free 
energy density with respect to the trial functions go{p), gi{p) and g{p) yields three saddle-point 
equations, 



9oip) 



(gip) - gi{p) + Xoigiip) - go{p))f 

gi{p) -gojp) 

{g{p) - gi{p)) {g{p) - giip) + Xo{gi{p) - go{p))) 



'^u [go] 

= 2m [(^1 - ^0)] 



(6.14) 
(6.15) 
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1 



9{p) -9i{p) 
The last equation can be rewritten as follows 



/ + r + (3(7 - 2u) [g] + 2u [gi 



(6.16) 



gip) -9i{p) 



p^ + Ai 



Ge(p; Ai) 



(6.17) 



where the unknown parameter Ai is defined by 



Ai = r + (3(7 - 2u) [GJ + 3(7 [gi] 
From eqs. (|6.15 ) and ( |6.14 ) one finds 



(6.18) 



9i{p) - 9o{p) = qiGc{p;Xi)Gc{p]{Xi-qiXo)) 

9o{p) = qoGl {p; (Ai - qiXo)) , 
and the unknown parameters qi and go are defined by 



(6.19) 
(6.20) 



[a(Ai)a(Ai - qiXo)] 



go (2m Gl{\i - qiXa 



- 1 



1 

2^ 
. 



(6.21) 



(6.22) 



One more saddle-point equation is obtained by taking the derivative of the free energy, eq.( |6.13|) , 
with respect to the parameter Xo', it reads 



[InG'c(Ai)] - [lnG'c(Ai - qiXo)] + qiXo [Gc(Ai - qiXo)] = —{qix 



4m 



(6.23) 



In this way, we have obtained four equations, ( |6.18D , (|6.21|) - (|6.23|) , for four unknown parameters, 
Ai, gi, go and Xo- 

Introducing the notation Ai — qiXo = A2 and taking into account the definition of the Green 
function Gc, eq.( |6.17|) , the equations ( |6.21|) and ( |6.23| ) can be represented as follows: 



[a(A2)]-[a(Ai)]--(Ai-A2) 



dX ([G,iX2)] - [^(A)] - i-(A - A2) 



(6.24) 



(6.25) 



One can easily check that the function V'(A) = [G'c(A2)] — [Gc(A)] — ^{X — A2) has a unique 
extremum at A = Aq (defined by [G^(Ao)] = 2^). Therefore, the two equations ( |6.24| ) and ( |6.25| ) can 
be simultaneously satisfied only when Ai = A2, which means that qiXo = 0. In either of the two 
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cases, gi = or Xo = 0, we come back to the replica-symmetric structure of the trial replica matrix 

Gab- 

Quite similar (although more cumbersome) calculations show that there are also no solutions 
with more than one step of RSB. Thus, within the present approximation we cannot have solutions 
with broken replica symmetry in the spin- glass phase. 



7 Appendix B 

The saddle-point solutions for non-ferromagnetic states (m = 0) are defined by two order param- 
eters A and /i. The stability of these states is determined by the signs of the eigenvalues of the 
corresponding Hessian: 



/ ^ 



\ 



92/ 



dfidX 



(7.1) 



As usual in the replica theory, one should take into account that in the process of taking the limit 
n — *• the minima of the replica free energy (at n > 1) can turn into maxima (at n < 1). The 
physically relevant states correspond of course to minima of the replica free energy before taking 
the limit n — 0. Around these states all the eigenvalues of the corresponding Hessian (composed 
of the second derivatives of the free energy with respect to all replica components of the order 
parameters) must be positive. However, if we consider the expression for the free energy where the 
limit n — >^ is already taken, the situation can change. In particular, within the replica-symmetric 
ansatz for the matrix fiab, eg. (|2.271) , the total number of off-diagonal elements (all equal to — /i) is 
equal to n{n — 1), and this number becomes negative for n < 1. Therefore, the physically relevant 
state defined by the free energy in eq.( p.34| ) (where the limit n — > is already taken) must be 
the maximum with respect to the parameter /i. On the other hand, the total number of diagonal 
elements of the matrix fiab is equal to n, and it remains positive in the limit n ^ 0. Therefore, the 
physically relevant extremum of the free energy, eq.( |2.34| ), must be the minimum with respect to 
the parameter A. Thus, in terms of the Hessian, eq. (|7. Ij) , a physically relevant saddle-point solution 
must correspond negative (corresponding to the parameter /i) and one positive (corresponding to 
the parameter A) eigenvalue. 

By taking the derivatives of the free energy, eq. (|2.34|) , with respect to and A and using the 
saddle-point equations (|3^) and ( p.7|) ), we find (for m = 0) 
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(7.3) 
(7.4) 
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For the paramagnetic state (/i = 0) we obtain 



dij? 
dfidX 
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[l + {3g-2u) [g']) 



The two eigenvalues of the Hessian, eg. (|7. 1| ) , are then 
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(7.6) 
(7.7) 
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(7. 



One of these eigenvalues is always positive, while the other one is negative when 2u [G"^] > 1; in this 
case the paramagnetic solution is stable. On the other hand, when 2u [G"^] < 1 the second eigenvalue 
is also positive, and the paramagnetic solution is unstable. Using eq.( p.2|) we find that the condition 
2u [G^] > 1 is equivalent to A(r) < Aq (eq.( p.lOD ), which is satisfied when r > Tgg. Thus, we can 
conclude that the paramagnetic solution is stable only when r > Tgg, while for r < Tgg it becomes 
unstable. 

We study now the stability of the spin-glass solution. Substituting eq.( p.8|) and [G^] = G/2Xo 
into eqs.(|7]|)-(|7j) leads to 
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The two corresponding eigenvalues are then 



(7.9) 
(7.10) 
(7.11) 



1 ± 



1 + 



3Gg 



r/i 



(7.12) 



where f'^x and f'^^ are short-hand notations for d'^f/dfi'^ and d^f/dX"^, respectively. From the above 
expression, one can conclude that in the temperature region r < r^g where the spin-glass solution 
is physical (/u(t) > 0, eq.( 3.11| )), there is always one positive and one negative eigenvalue, which 
means that the spin-glass solution is stable. 

Thus for T > Tgg the only stable state of the system is paramagnetic, while for r < Tsg the 
paramagnetic solution becomes unstable and the stable state of the system is the spin-glass phase. 
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8 Appendix C 



For the spin-glass state in dimensions D = 4 — e (e <^ 1) one again obtains the solutions, 
eqs(|378| ) , (|3.9| ), where instead of eqs. (|3.1| ),( ^^ one now has 



p\<l {2Tl)Dp^ + \ \ t ^ 



C 



A<1 



Then, instead of eg. ( |3.1CI| ), eg. (pl8|) leads to 



A. 



4Cn(l-f) 
e + ACu 



•2) 



In the limit e -C m <C 1 this result reduces back to eg. (|3.1CI| ), so that one eventually recovers all the 
solutions corresponding to the case Z) = 4 described in section 3. 



Let us consider the other limit 



In this case eg. ( ^I3|) yields 



M < e < 1 
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'AC 



-u 



\ e 



15) 



Thus, one can easily see that when u <^ e <^ 1 one recovers all the solutions described in section 3, 



in which the value of Aq is now given by eg.( |8.5| ). In particular, it is this value of Aq that controls 
the temperature range of the spin-glass phase separating the paramagnetic and the ferromagnetic 
phases. 

As in section 5.1 where we discussed the validity of the Gaussian variational approach in di- 
mensions D = 4, one should study in place of of egs.( |5.2| ),( ^.3| ), the following RG eguations: 



d_ 



giO = eg-6C {^g - Au) g + O ig'; g\ gu^; 



d_ 



u{^) = eu-AC {3g - 4u) u + O ( ^f^; g^u; gu^] 



17) 



These eguations show that the higher-order corrections of the perturbation theory remain irrelevant 
at spatial scales R bounded by the condition 



6C|3^o -4mo|^' < e 
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which corresponds to the temperature scale 



T > T* ~ I I . (8.9) 

Thus, as done in section 5.1, to guarantee that the higher-order corrections do not destroy the 
approximation made in the present approach, one needs Aq ^ r*, or 

One can easily verify that this condition is satisfied when 

H < * < 11 (8.11) 

9 Mo 9 ^ ^ 

and e -C 1. In other words, we find the same restriction on the parameters g and u as in D = 4, 
but in addition one needs the value oi e — 4 — D to he small. 
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